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Abstract 

We review the formalism of the effective average action in quantum field theory which 
corresponds to a coarse grained free energy in statistical mechanics. The associated ex- 
act renormalization group equation and possible nonperturbative approximations for its 
solution are discussed. This is applied to QCD where one observes the consecutive emer- 
gence of mesonic bound states and spontaneous chiral symmetry breaking as the coarse 
graining scale is lowered. We finally present a study of the phenomenological importance 
of non-renormalizable terms in the effective linear meson model. 

1 Effective average action 

Quantum chromodynamics (QCD) describes qualitatively different physics at different length 
scales. At short distances the relevant degrees of freedom are quarks and gluons which can 
be treated perturbatively. At long distances we observe hadrons and an essential part of the 
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dynamics can be encoded in the masses and interactions of mesons. Any attempt to deal with 
this situation analytically and to predict the meson properties from the short distance physics 
(as functions of the strong gauge coupling and the current quark masses m^) has to bridge 
the gap between two qualitatively different effective descriptions. Two basic problems have to 
be mastered for an extrapolation from short distance QCD to mesonic length scales: 

• The effective couplings change with scale. This does not only concern the running gauge 
coupling, but also the coefficients of non-renormalizable operators as, for example, four 
quark operators. Typically, these non-renormalizable terms become important in the 
momentum range where ag is strong and deviate substantially from their perturbative 
values. Consider the four-point function which obtains after integrating out the gluons. 
For heavy quarks it contains the information about the shape of the heavy quark potential 
whereas for light quarks the complicated spectrum of light mesons and chiral symmetry 
breaking are encoded in it. At distance scales around Ifm one expects that the effective 
action resembles very little the form of the classical QCD action which is relevant at short 
distances. 

• Not only the couplings, but even the relevant variables or degrees of freedom are different 
for long distance and short distance QCD. It seems forbiddingly difficult to describe 
the low-energy scattering of two mesons in a language of quarks and gluons only. An 
appropriate analytical field theoretical method should be capable of introducing field 
variables for composite objects such as mesons. 

A conceptually very appealing idea for our task is the block-spin action 0. It realizes 
that physics with a given characteristic length scale / is conveniently described by a functional 
integral with an ultraviolet (UV) cutoff A for the momenta. Here A should be larger than 
but not necessarily by a large factor. The Wilsonian effective action S]^ replaces then the 
classical action in the functional integral. It obtains by integrating out the fluctuations with 
momenta ^ A^. An exact renormalization group equation describes how changes 

with the UV cutoff A. 

We will use here the somewhat different but related concept of the effective average action 
Ffc which, in the language of statistical physics, is a coarse grained free energy with coarse 
graining scale k. The effective average action is based on the quantum field theoretical concept 
of the effective action |^ F which obtains by integrating out all quantum fluctuations. The 
effective action contains all information about masses, couplings, form factors and so on, since it 
is the generating functional of the IPI Green functions. The field equations derived from F are 
exact including all quantum effects. For a field theoretical description of thermal equilibrium 
this concept is easily generalized to a temperature dependent effective action which includes now 
also the thermal fluctuations. In statistical physics F describes the free energy as a functional 
of some convenient (space dependent) order parameter, for instance the magnetization. In 
particular, the behavior of F for a constant order parameter (the effective potential) specifies 
the equation of state. The effective average action F^ is a simple generalization of the effective 
action, with the distinction that only quantum fluctuations with momenta ^ k"^ are included. 
This can be achieved by introducing in the functional integral defining the partition function (or 
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the generating functional for the n-point functions) an exphcit infrared cutoff ~ k. Typically, 
this IR-cutoff is quadratic in the fields and modifies the inverse propagator, for example by 
adding a mass-like term ~ k^. The effective average action can then be defined in complete 

analogy to the effective action (via a Legendre transformation of the logarithm of the partition 
function). The mass-like term in the propagator suppresses the contributions from the small 
momentum modes with ^ k"^ and accounts effectively only for the fluctuations with 
^ k\ 

Following the behavior of F^ for different k is hke looking at the world through a microscope 
with variable resolution: For large k one has a very precise resolution ~ k~^ but one also studies 
effectively only a small volume ~ k'^. Taking in QCD the coarse graining scale k much larger 
than the confinement scale guarantees that the complicated nonperturbative physics does not 
play a role yet. In this case, F^ will look similar to the classical action, typically with a running 
gauge coupling cvahiated at the scale k. (This does not hold for Green functions with much 
larger momenta S> fc^ where the relevant IR cutoff is p and the effective coupling as{p)-) 
For lower k the resolution is smeared out and the detailed information of the short distance 
physics can be lost. (Again, this does not concern Green functions at high momenta.) On the 
other hand, the "observable volume" is increased and long distance aspects such as collective 
phenomena become now visible. In a theory with a physical UV cutoff A we may associate 
Fa with the classical action S since no fluctuations are effectively included. By definition, the 
effective average action equals the effective action for A: = 0, Fq = F, since the infrared cutoff 
is absent. Thus F^ interpolates between the classical action S and the effective action F as /c is 
lowered from A to zero. The ability to follow the evolution to A; —> is equivalent to the abihty 
to solve the quantum field theory. 

For a formal description we will consider in the first two sections a model with real scalar 
fields x"" (the index a labeling internal degrees of freedom) in d Euclidean dimensions with 
classical action 5'[x]. We define the generating functional for the connected Green functions by 

Wk[J] ^InJ Pxexp [Skbc] + J o?'^xJ„(x)x"(x)} (1) 
where we have added to the classical action an IR cutoff A^S 

S,[x] = S[x] + AS,[x] (2) 
which is quadratic in the fields and best formulated in momentum space 

A-S.[x] = IJ ^,Rk{q')Xa{-q)x''{q) ■ (3) 

Here are the usual scalar sources introduced to define generating functionals and Rk{q^) 
denotes an appropriately chosen (see below) IR cutoff function. We require that Rk{q^) vanishes 
rapidly for q^ ^ k"^ whereas for q^ <^ k"^ it behaves as Rk{q'^) — k"^ . This means that all Fourier 
components of x" with momenta smaller than the IR cutoff k acquire an effective mass mefr — k 
and therefore decouple while the high momentum components of x" should not be affected by 
Rk- The classical fields 

. ^ ^ (4) 
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depend now on k. In terms of Wk the effective average action is defined via a Legendre transform 

= -Wk[J] + J d^xUx)<!>-{x) - AS^m . (5) 

In order to define a reasonable coarse grained free energy we have subtracted in the infrared 
cutoff piece. This guarantees that the only difference between and F is the effective IR cutoff 
in the fluctuations. Furthermore, this has the consequence that F^ does not need to be convex 
whereas a pure Legendre transform is always convex by definition. (The coarse grained free 
energy becomes convex only for k 0.) This is very important for the description of phase 
transitions, in particular first order ones. One notes 

limi?fc(g2) = limFfc[<l>] = F[$] 

lim i?fc(g2) = oo ^ limF,[<I>] = 5[<l>] 
for a convenient choice of Rk like 

Rk{q^) = Zkq^—p^^^—^^^^ . (7) 

Here denotes the wave function renormalization to be specified below and we will often use 
for Rk the limit A ^ oo 

= iMzr^ ■ (8) 

We note that the property Fa = 5 is not essential since the short distance laws may be 
parameterized by Fa as well as by S. In addition, for momentum scales much smaller than A 
universality implies that the precise form of Fa is irrelevant, up to the values of a few relevant 
renormalized couplings. 

A few properties of the effective average action are worth mentioning: 

1. All symmetries of the model that are respected by the IR cutoff A^S" are automatically 
symmetries of F^. In particular, this concerns translation and rotation invariance and one 
is not plagued by many of the problems encountered by a formulation of the block-spin 
action on a lattice. 

2. In consequence, F^, can be expanded in terms of invariants with respect to these symme- 
tries with couplings depending on k. For the example of a scalar theory one may use a 
derivative expansion (p = $"^$0/2) 

Ffe = 1 d'^x |f/fc(p) + ]^Zk{p)d^^ad^^'' + . . .} (9) 

and expand further in powers of p 

1_ 2 1_ 3 

Uk{p) = -^hip-poik)) + g7fc(p-Po(^)) +•■■ 

Zk{p) = Zk{po{k)) + Z',{po{k)){p-po{k)) + .... (10) 

We see that F^ describes infinitely many running couplings. {Zk in (^ can be identified 
with Zk{po).) 
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3. There is no problem incorporating chiral fermions since a chirally invariant cutoff Rk can 
be formulated JTO] . 



Gauge theories can be formulated along similar lines [|TT]-[|I6| even though A^S may not 
be gauge invariant. In this case the usual Ward identities receive corrections for which 
one can derive closed expressions These corrections vanish for k ^ 0. 



The high momentum modes are very effectively integrated out because of the exponential 
decay of Rk for ^ fc^. Nevertheless, it is sometimes technically easier to use a cutoff 
without this fast decay property (e.g. Rk ~ k"^ or Rk ~ k^/q^). In the latter case one has 
to be careful with possible remnants of an incomplete integration of the short distance 
modes. Also our cutoff does not introduce any non-analytical behavior as would be the 
case for a sharp cutoff [0. 

Despite a similar spirit and many analogies there remains also a conceptual difference 
to the Wilsonian effective action . The Wilsonian effective action describes a set 
of different actions (parameterized by A) for one and the same theory — the n-point 
functions are independent of A and have to be computed from by further functional 
integration. In contrast, Tk describes the effective action for different theories — for any 
value of k the effective average action is related to the generating functional of a theory 
with a different action Sk = S + AkS. The n-point functions depend on k. The Wilsonian 
effective action does not generate the IP/ Green functions [1^ . 



7. Because of the incorporation of an infrared cutoff, Tk is closely related to an effective 
action for averages of fields [0, where the average is taken over a volume ~ k*^. 



2 Exact renormalization group equation 



The dependence of Tk on the coarse graining scale k is governed by an exact renormalization 
group equation (ERGE) pf 



a,r,[$] = -Tr^ T^m+Rk 



dtRk 



Here t = ln(fc/A) with some arbitrary momentum scale A, and the trace includes a momentum 



integration as well as a summation over internal indices, Tr = / -^i^^ Ha- The second functional 
derivative F^^-* denotes the exact inverse propagator 

5'T, 



T 



(2) 



(12) 



The flow equation (|TTD can be derived from (^ in a straightforward way using 

dtTkU = -dtWk\j-dAkSm 
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^Tr{a,wf} (13) 



and 



6Ja{-q)6Mq') 6M-q')6^,{q") 



SaJqq" ■ (14) 



It has the form of a renormahzation group improved one-loop expression 0. Indeed, the 
one-loop formula for reads 

T,m = Sm + ^TT\n{S^'\^ + Rk) (15) 

with 5**^^-* the second functional derivative of the classical action, similar to (|13). (Remember 
that S**^^-* is the field dependent classical inverse propagator. Its first and second derivative with 
respect to the fields describe the classical three- and four-point vertices, respectively.) Taking 
a t-derivative of eq. (0) gives a one-loop flow equation very similar to (|1T]) with T^^^ replaced 
by S''-^-'. It may seem surprising, but it is nevertheless true, that the renormahzation group 
improvement S*^^-* T^^^ promotes the one-loop flow equation to an exact nonperturbative 
flow equation which includes the effects from all loops as well as all contributions which are 
non-analytical in the couplings like instantons, etc.! For practical computations it is actually 
often quite convenient to write the flow equation (pT]) as a formal derivative of a renormahzation 
group improved one-loop expression 

aiFfc = iTr5tln(Fl')+/2,) (16) 

with dt acting only on and not on T^^\ i.e. dt = (dRk/dt) (d/dRk). Flow equations for 
n-point functions follow from appropriate functional derivatives of (|TT|) or (|TBp with respect to 
the fields. For their derivation it is sufficient to evaluate the corresponding one-loop expressions 
(with the vertices and propagators derived from F^) and then to take a formal 9f-derivative. (If 
the one-loop expression is finite or properly regularized the 9j-derivative can be taken after the 
evaluation of the trace.) This permits the use of (one-loop) Feynman diagrams and standard 
perturbative techniques in many circumstances. Most importantly, it establishes a very direct 
connection between the solution of flow-equations and perturbation theory. If one uses on the 
right hand side of eq. ([TT|) a truncation for which the propagator and vertices appearing in 



are replaced by the ones derived from the classical action, but with running /c-dependent 
couplings, and then expands the result to lowest non-trivial order in the coupling constants one 
recovers standard renormahzation group improved one-loop perturbation theory. The formal 
solution of the flow equation can also be employed for the development of a systematically 
resummed perturbation theory 
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For a choice of the cutoff function similar to (|]) the momentum integral contained in the 
trace on the right hand side of the flow equation is both infrared and ultraviolet finite. In- 
frared finiteness arises through the presence of the infrared regulator ~ R^. We note that all 
eigenvalues of the matrix + -Rfc must be positive semi-definite. The proof follows from the 
observation that the functional F^ + A^S" is convex since it is obtained from by a Legen- 
dre transform. On the other hand, ultraviolet finiteness is related to the fast decay of dtRk 
for ^ k^. This expresses the fact that only a narrow range of fluctuations with ~ k'^ 
contributes effectively if the infrared cutoff k is lowered by a small amount Q Since the flow 
equation is manifestly flnite this can be used to deflne a regularization scheme. The "ERGE- 
scheme" is specifled by the flow equation, the choice of Rk and the "initial condition" Fa. 
This is particularly important for gauge theories where other regularizations in four dimensions 
and in the presence of chiral fermions are difficult to construct [|. We note that in contrast 
to previous versions of exact renormalization group equations there is no need in the present 
formulation to construct an ultraviolet momentum cutoff — a task known to be very difficult 
in non-Abelian gauge theories. 

Despite the conceptual differences between the Wilsonian effective action and the ef- 
fective average action F^ the exact flow equations describing the A-dependence of S]^ and the 
^-dependence of F^ are simply related. Polchinski's continuum version of the Wilsonian flow 
equation can be transformed into (0) by means of a Legendre transform and a suitable 
variable redeflnition | |2(J[ | . 

Even though intuitively simple, the replacement of the (RG-improved) classical propagator 
by the full propagator turns the solution of the flow equation (|lT]) into a difficult mathematical 
problem: The evolution equation is a functional differential equation. Once F^ is expanded 
in terms of invariants (e.g. (p!oD) this is equivalent to a coupled system of non-linear 
partial differential equations for infinitely many couplings. General methods for the solution 
of functional differential equations are not developed very far. They are mainly restricted to 
iterative procedures that can be applied once some small expansion parameter is identified. 
This covers usual perturbation theory in the case of a small coupling, the 1/A^-expansion or 
expansions in the dimensionality A — d or 2 — d. It may also be extended to less familiar 
expansions like a derivative expansion which is related in critical three dimensional scalar 
theories to a small anomalous dimension rj. In the absence of a clearly identified small parameter 
one nevertheless needs to truncate the most general form of F^ in order to reduce the infinite 
system of coupled differential equations to a (numerically) manageable size. This truncation is 
crucial. It is at this level that approximations have to be made and, as for all nonperturbative 
analytical methods, they are often not easy to control. The challenge for nonperturbative 
systems like low momentum QCD is to find fiow equations which (a) incorporate all the relevant 
dynamics such that neglected effects make only small changes, and (b) remain of manageable 
size. The difficulty with the first task is a reliable estimate of the error. For the second task 

^If for some other choice of Rk the right hand side of the flow equation would not remain UV finite this 
would indicate that the high momentum modes have not yet been integrated out completely in the computation 
of Tfc. 

^For gauge theories Fa has to obey appropriately modified Ward identities. In the context of perturbation 
theory a first proposal how to regularize gauge theories by use of flow equations can be found in |12, 03) . 
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the main limitation is a practical restriction for numerical solutions of differential equations to 
functions depending only on a small number of variables. The existence of an exact functional 
differential flow equation is a very useful starting point and guide for this task. At this point 
the precise form of the exact flow equation is quite important. Furthermore, it can be used for 
systematic expansions through enlargement of the truncation and for an error estimate in this 
way. Nevertheless, this is not all. Usually, physical insight into a model is necessary to device 
a useful nonperturbative truncation! 

So far, two complementary approaches to nonperturbative truncations have been explored: 
an expansion of the effective Lagrangian in powers of derivatives (p = 

Tf^m = J d^x \Uk{p) + + h^k{p)d^pd>'p + (17) 

or one in powers of the fields 



oo 




If one chooses $o as the fc-dependent VEV of $, the series (O) starts effectively at n = 2. 

in] 

The flow equations for the IP I n-point functions V\ are obtained by functional differentiation 
of eq. (pUD 0. The formation of mesonic bound states, which typically appear as poles in the 
(Minkowskian) four quark Green function, is most efficiently described by expansions like (p!^). 
This is also the form needed to compute the nonperturbative momentum dependence of the 
gluon propagator and the heavy quark potential |T5|, |TB[. On the other hand, a parameterization 
of Ffc as in (|1^ seems particularly suited for the study of phase transitions. The evolution 
equation for the average potential Uk follows by evaluating (|l^ for constant $. In the limit 
where the ^-dependence of is neglected and = one finds for the 0(A^)-symmetric 
scalar model 

an ( \ = - [ dRk ( N-l 1 \ 

' "^P^ 2J {2'nY dt \Z,q^ + R, + U', Z,q^ + R,{q) + U', + 2pU';) ^ ^ 

with t/^ = etc. One observes the appearance of p-dependent mass terms in the effective 
propagators of the right hand side of eq. ([19|) . Once rj^ = —dtlnZj, is determined in terms 
of the couplings parameterizing Uk this is a partial differential equation Q for a function Uk 
depending on two variables k and p which can be solved numerically |2^-||2^. (The Wilson- 
Fisher fixed point relevant for a second order phase transition [d = 3) corresponds to a scaling 
solution ||2^, ^ where dtUk = 0.) A suitable truncation of a flow equation of the type (p!7|) will 



play a central role in the description of chiral symmetry breaking below. 



■^Such flow equations have been discussed earlier from a somewhat different viewpoint They can also be 
interpreted as a differential form of Schwinger-Dyson equations 

*In the sharp cutoff limit and for rjk = the flow equation for Uk coincides with the one corresponding to the 
Wilsonian effective action Q . Because of non-analyticities in the kinetic terms [Q an analysis of the derivative 
terms is difficult for a sharp cutoff. 
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It should be mentioned at this point that the weakest point in the ERGE approach seems 
to be a rehable estimate of the truncation error in a nonperturbative context. This problem is 
common to all known analytical approaches to nonperturbative phenomena and appears often 
even within systematic (perturbative) expansions. One may hope that the existence of an exact 
flow equation could also be of some help for error estimates. An obvious possibility to test a 
given truncation is its enlargement to include more variables — for example, going one step 
higher in a derivative expansion. This is similar to computing higher orders in perturbation 
theory and limited by practical considerations. As an alternative, one may employ different 
truncations of comparable size — for instance, by using different definitions of retained cou- 
plings. A comparison of the results can give a reasonable picture of the uncertainty if the used 
set of truncations is wide enough. In this context we should also note the dependence of the 
results on the choice of the cutoff function Rk{q)- Of course, for — the physics should 
not depend on a particular choice of Rk and, in fact, it does not for full solutions of (|Tl|). 
Different choices of Rk just correspond to different trajectories in the space of effective actions 
along which the unique IR limit r[$] is reached. Once approximations are used to solve the 
ERGE (|TT|), however, not only the trajectory but also its end point will depend on the precise 
definition of the function Rk- This is very similar to the renormalization scheme dependence 
usually encountered in perturbative computations of Green functions. One may use this scheme 
dependence as a tool to study the robustness of a given approximation scheme. 

Before applying a new nonperturbative method to a complicated theory like QCD it should 
be tested for simpler models. A good criterion for the capability of the ERGE to deal with 
nonperturbative phenomena concerns the critical behavior in three dimensional scalar theories. 
In a first step the well known results of other methods for the critical exponents have been 
reproduced within a few percent accuracy pGf. The ability of the method to produce new 



results has been demonstrated by the computation of the critical equation of state for Ising 

This has been 
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or for the Abelian Higgs model relevant 



and Heisenberg models p3| which has been verified by lattice simulations 
extended to first order transitions in matrix models 

for superconductors p9|, ^ . Analytical investigations of the high temperature phase transitions 
in d = 4 scalar theories (O(A^)-models) have correctly described the second order nature of the 
transition 
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in contrast to earlier attempts within high temperature perturbation theory. 
For an extension of the flow equations to Abelian and non-Abelian gauge theories we refer 
the reader to |Tl|-|T^. The other necessary piece for a description of low-energy QCD, namely 
the transition from fundamental (quark and gluon) degrees of freedom to composite (meson) 

We will describe the most 



fields within the framework of the ERGE can be found in 
important aspects of this formalism for mesons below. 
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3 Chiral symmetry breaking in QCD 

The strong interaction dynamics of quarks and gluons at short distances or high energies is 
successfully described by quantum chromo dynamics (QCD). One of its most striking features 



is asymptotic freedom which makes perturbative calculations reliable in the high energy 



regime. On the other hand, at scales around a few hundred MeV confinement sets in. As 
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a consequence, the low-energy degrees of freedom in strong interaction physics are mesons, 
baryons and glueballs rather than quarks and gluons. When constructing effective models for 
these IR degrees of freedom one usually relies on the symmetries of QCD as a guiding principle, 
since a direct derivation of such models from QCD is still missing. The most important sym- 
metry of QCD, its local color SU{3) invariance, is of not much help here, since the IR spectrum 
appears to be color neutral. When dealing with bound states involving heavy quarks the so 
called "heavy quark symmetry" may be invoked to obtain approximate symmetry relations 



between IR observables We will rather focus here on the light scalar and pseudoscalar 

meson spectrum and therefore consider QCD with only the light quark flavors u, d and s. To 
a good approximation the masses of these three flavors can be considered as small in com- 
parison with other typical strong interaction scales. One may therefore consider the chiral 
limit of QCD (vanishing current quark masses) in which the classical QCD Lagrangian does 
not couple left- and right-handed quarks. It therefore exhibits a global chiral invariance under 
Ul{N) X Ur{N) = SUiiN) X SUr{N) x f/y (1) x f/^(l) where denotes the number of massless 
quarks (A^ = 2 or 3) which transform as 

^i^i^^V- ^LipL-, ULeUL{N). (20) 

Even for vanishing quark masses only the diagonal SUv{N) vector-like subgroup can be ob- 
served in the hadron spectrum ("eightfold way"). The symmetry SUl{N) x SUr{N) must 
therefore be spontaneously broken to SUv{N) 

SUl{N) X SUr{N) ^ SUlMN) = SUv{N) . (21) 

Chiral symmetry breaking is one of the most prominent features of strong interaction dynam- 



ics and phenomenologically well established |3J], even though a rigorous derivation of this 



phenomenon starting from first principles is still missing. In particular, the chiral symmetry 
breaking (|2T|) predicts for = 3 the existence of eight light parity-odd (pseudo-)Goldstone 
bosons: 7r°, vr^, K^, K , K"^ and Tj. Their comparably small masses are a consequence of the 
explicit chiral symmetry breaking due to small but non-vanishing current quark masses. The 
axial Abelian subgroup Ua{^) = Ul^r{1) is broken in the quantum theory by an anomaly of the 
axial-vector current. This breaking proceeds without the occurrence of a Goldstone boson . 
Finally, the Uv{^) = Ul+r{^) subgroup corresponds to baryon number conservation. 

The light pseudoscalar and scalar mesons are thought of as color neutral quark-antiquark 
bound states ^'^'^ ~ ^l'^r? C'^b = 1, . . . ,N, which therefore transform under chiral rotations 
jQl) as 

$ . (22) 



Hence, the chiral symmetry breaking pattern (^I]) is realized if the meson potential develops a 
VEV 



= ao5'^^ ; ao^O. (23) 
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One of the most crucial and yet unsolved problems of strong interaction dynamics is to derive 
an effective field theory for the mesonic degrees of freedom directly from QCD which exhibits 
this behavior. 



4 A semi— quantitative picture 

Before turning to a quantitative description of chiral symmetry breaking using flow equations 
it is perhaps helpful to give a brief overview of the relevant scales which appear in relation to 
this phenomenon and the physical degrees of freedom associated to them. Some of this will 
be explained in more detail in the remainder of this talk whereas other parts are rather well 
established features of strong interaction physics. 

At scales above approximately 1.5 GeV, the relevant degrees of freedom of strong interactions 
are quarks and gluons and their dynamics appears to be well described by perturbative QCD. 
At somewhat lower energies this changes dramatically. Quark and gluon bound states form and 
confinement sets in. Concentrating on the physics of scalar and pseudoscalar mesons Q there 
are three important momentum scales which appear to be rather well separated: 

• The compositeness scale at which mesonic ipip bound states form due to the increasing 



strength of the strong interaction. It will turn out to be somewhere in the range (600 — 
700) MeV. 



• The chiral symmetry breaking scale k^sB at which the chiral condensate y}) tp'^J or 
will assume a no n- vanishing value, therefore breaking chiral symmetry according to (pl|). 
This scale is found to be around (400 — 500) MeV. Below it the quarks acquire constituent 
masses Mg ~ 350 MeV due to their Yukawa coupling to the chiral condensate (^3]). 

• The confinement scale Aqcd — 200 MeV which corresponds to the Landau pole in the 
perturbative evolution of the strong coupling constant a^. In our context, this is the scale 
where deviations of the effective quark propagator from its classical form and multi-quark 
interactions not included in the meson physics become of crucial importance. 

For scales k in the range k^sB ^ k ^ kq> the most relevant degrees of freedom are mesons and 
quarks. Typically, the dynamics in this range is dominated by the strong Yukawa coupling 
h between quarks and mesons: h"^ {k) / {Att) ^ as{k)- One may therefore assume that the 
dominant QCD effects are included in the meson physics and consider a simple model of quarks 
and mesons only. As one evolves to scales below k^sB the Yukawa coupling decreases whereas as 
increases. Of course, getting closer to Aqcd it is no longer justified to neglect the QCD effects 
which go beyond the dynamics of effective meson degrees of freedom Q. On the other hand, 
the final IR value of the Yukawa coupling h is fixed by the typical values of constituent quark 
masses Mq ~ 350 MeV to be h'^/{ATT) ~ 4.5. One may therefore speculate that the domination 

^One may assume that all other bound states are integrated out. We will comment on this issue below, 
formalism for the inclusion of "residual" gluon contributions to the flow equation in a situation where 
part of the gluon effects are already described by mesonic bound states can be found in p5| . 
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of the Yukawa interaction persists down to scales k ~ Mq at which the quarks decouple from 
the evolution of the mesonic degrees of freedom altogether due to their mass. Of course, details 
of the gluonic interactions are expected to be crucial for an understanding of quark and gluon 
confinement. Strong interaction effects may dramatically change the momentum dependence 
of the quark and gluon propagators for k around Aqcd- Yet, as long as one is only interested 
in the dynamics of the mesons one is led to expect that these effects are quantitatively no too 
important. Because of the effective decoupling of the quarks and therefore the whole colored 
sector the details of confinement have only little influence on the mesonic flow equations for 
^ ~ Aqcd- 

Let us imagine that we integrate out the gluon degrees of freedom while keeping an effective 
infrared cutoff kp ~ 1.5 GeV in the quark propagators. The exact flow equation to be used 
for this purpose obtains by keeping for k < kp the infrared cutoff R for the quarks fixed 
while lowering the one for the gluons to zero 0. Subsequently, the gluons are eliminated by 
solving the field equations for the gluon fields as functionals of the quarks. This will result 
in a non-trivial momentum dependence of the quark propagator and effective non-local four 
and higher quark interactions. Because of the infrared cutoff the resulting effective action for 
the quarks resembles closely the one for heavy quarks (at least for Euclidean momenta). The 
dominant effect is the appearance of an effective quark potential (similar to the one for the 
charm quark) which describes the effective four quark interactions. For the effective quark 
action at kp we only retain this four quark interaction in addition to the classical two-point 
function, while neglecting ra-point functions involving six and more quarks. Details of the quark 
two-point function for ^ ^qcb will important as long as A; > Aqcd, and similarly 

for the four-point function. For k < Aqcd the quarks decouple from the mesonic sector as 
discussed above. For typical momenta larger than Aqcd a reliable computation of the effective 
quark action should be possible by using in the quark-gluon flow equation relatively simple 
truncations |l^, [16[ . 

We next have to remove the infrared cutoff for the quarks, k ^ 0. This task can be carried 
out by means of the exact flow equation for quarks only, starting at kp with an initial value 
Tfepi'?/'] as obtained after integrating out the gluons. For fermions the trace in eq. (0) has 
to be replaced by a supertrace in order to account for the minus sign related to Grassmann 
variables |T^. A first investigation in this direction has used a truncation with a chirally 
invariant four quark interaction whose most general momentum dependence was retained 



+ 



X 



d'^p —i 



(27r)4 

^''k'\pi^P2,P3,PA) {[^'ai-pM{P2)] [^i{PA)^]{-P3 



n ) (27r)^ (5(pi +P2-P3- Pa) 



V'l(-Pl)75^»^(P2)] [V^6(P4)75^J(-P3)]} • 



(24) 



""One may also only partly integrate out the gluons by stopping this evolution at a scale somewhat above 
Aqcd- The remaining residual gluon fluctuation can then be included later following Jist . 
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Here i, j run from one to Nc which is the number of quark colors. The indices a, b denote different 
hght quark flavors and run from 1 to A^. The matrices m and m are hermitian and m + 77175 
forms therefore the most general quark mass matrix 0. The ansatz (p^ does not correspond 
to the most general chirally invariant four quark interaction. It neglects similar interactions in 
the p-meson and pomeron channels which are also obtained from a Fierz transformation of the 
heavy quark potential [0. With Vlq"^) the heavy quark potential in a Fourier representation, 



the initial value at kp = 1.5 GeV is given by {Z^^k = Z^^kip = —kp)) 

>^^lt\puP2,P3,P4)Z^X = l^iiP^-Psf) = ( ^^"%2 + ^o^2 • (25) 

For simplicity, the effective heavy quark potential is approximated here by a one gluon exchange 
term ~ as{kp) and a string tension A ~ 0.18 GeV^, but one may use as well a more precise 
form of the potential as determined from the phenomenology of charmonium or as computed 



from the flow equations with gluons |15, 16|. In the simplified ansatz (25) the string tension 
introduces a second scale in addition to kp and it becomes clear that the incorporation of gluon 
fluctuations is a crucial ingredient for the emergence of mesonic bound states. For a more 
realistic treatment of the heavy quark potential this scale is set by the running of or Aqcd- 
The evolution equation for the function A^'^'' can be derived from the fermionic version of 
eq. (p!T|) and the truncation p^). Since A^'^'' depends on six independent momentum invariants 
it is a partial differential equation for a function depending on seven variables and has to be 
solved numerically ||31|. The ansatz p5| ) corresponds to the t-channel exchange of a colored 
gluonic state and it is by far not clear that the evolution of A^'^'' will lead at lower scales to a 
momentum dependence representing the exchange of colorless mesonic bound states. Yet, at 
the compositeness scale 

~ 630 MeV (26) 

one finds an approximate factorization 

a2^(Pi,P2,P3,P4) = 9{pi,P2)G{s)g{p3,P4) + ■■■ (27) 

which indicates the formation of mesonic bound states. Here g{pi,P2) denotes the amputated 
Bethe-Salpeter wave function and G{s) is the mesonic bound state propagator displaying a 
pole-like structure in the s-channel if it is continued to negative s = {pi + P2)^ . The dots 
indicate the part of A^'^'' which does not factorize and which will be neglected in the following. 
In the limit where the momentum dependence of g and G is neglected we recover the four quark 
interaction of the Nambu-Jona-Lasinio model [^, It is therefore not surprising that our 
description of the dynamics for k < k^ will parallel certain aspects of the investigations of this 
model, even though we are not bound to the approximations used typically in such studies 
(large- A'^c expansion, perturbative renormalization group, etc.). 

It is clear that for scales ^ fc$ a description of strong interaction physics in terms of quark 
fields alone would be rather inefficient. Finding physically reasonable truncations of the effective 

*Our chiral conventions [ p^ where the hermitean part of the mass matrix is multiphed by 75 may be somewhat 
unusual but they are quite convenient for Euchdean calculations. 
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average action should be much easier once composite fields for the mesons are introduced. The 
exact renormalization group equation can indeed be supplemented by an exact formalism for 



the introduction of composite field variables or, more generally, a change of variables [31]. For 
our purpose, this amounts in practice to inserting at the scale kq, the identities ^ 



const 



Jva, 

X exp tr (4 - K],G -m\- O^G) 1 (a^ - GKa - ttia - GO) } 
1 = const / V(Th (2^ 



X exp ^ - i tr (alj - kIjG - m^H - C^^^^G) i (an - GKh - nin - GO^^^ 
1-2 G 

into the functional integral which formally defines the quark effective average action. Here Ka,h 
are sources for the collective fields (Ja,h which correspond in turn to the anti-hermitian and 
hermitian parts of the meson field They are associated to the fermion bilinear operators 
O ['?/'], whose Fourier components read 

o^'\{q) = - j -^,g{-P,P + q^\p)i,Mp + q)- (29) 

The choice of g{—p,p + q) as the bound state wave function renormalization and of G{q) as its 
propagator guarantees that the four-quark interaction contained in (^) cancels the dominant 
factorizing part of the QCD-induced non-local four-quark interaction (^4|), (pTj) . In addition, 
one may choose 



m 



Aab 



mab{0)g~\0,0)Z^,kM (30) 



such that the explicit quark mass term cancels out for q = 0. The remaining quark bilinear, 
which is proportional to m(g) — m{0)Z^ k^{0)g{—q,q)/[Z^ k^{q)g{0,0)] and breaks chiral sym- 
metry, vanishes for zero momentum and will be neglected in the following. Without loss of 
generality we can take m real and diagonal and m = 0. In consequence, we have replaced at 
the scale kq> the effective quark action (^4|) with ( pT]) by an effective quark meson action given 

by 

Tk = rfc-^|dStr($tj + jt$) 

Tk = I d^xUu{^,^^) (31) 



9 We use here the shorthand notation A'^GB = J -^Al{q)G''^{q)Bb{q). 

^°The fields aA,H associated to Ka.h by a Legendre transformation obey ct^ = — |(<i> — $^), — i(<i> + 
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At the scale /c$ the inverse scalar propagator is related to G{q) in eq. (^) by 

G-\q^) = 2m\k^) + 2Z$,fcJg)g2 . (32) 

This fixes the term in U^^ which is quadratic in $ to be positive, Uk^ = m^tr$^$ + . . .. The 
higher order terms in Uk^ cannot be determined in the approximation (^) since they correspond 
to terms involving six or more quark fields. The initial value of the Yukawa coupling corresponds 
to the "quark wave function in the meson" in eq. (pT)), i.e. 

hk^-q.q-p) = g{-q,q-p) (33) 

which can be normalized with /ifc,j,(0,0) = (7(0,0) = 1. We observe that the explicit chiral 
symmetry breaking from non-vanishing current quark masses appears now in the form of a 
meson source term with 

J = 2m^Z^^fc^(0)5f^^(0,0) {mab + irhab) = 2Z^^k^rrfdmg{mu, rrid, rris) . (34) 

This induces a non-vanishing ($) and an effective quark mass Mg through the Yukawa coupling. 
Spontaneous chiral symmetry breaking can be described in this language by a non-vanishing 
($) in the limit j — > 0. We note that the current quark mass nig and the constituent quark mass 
Mq ~ h ($) are identical only for a quadratic potential ~ tr Because of spontaneous chiral 
symmetry breaking the constituent quark mass Mg can differ from zero even for rrig = 0. One 
may check that by solving the field equation for $ as a functional of ip, ip (with Uk = WT^ tr $^$) 
one recovers from (^) the effective quark action (pi)). Note that in (|3l|) is chirally invariant. 
The flow equation for Tk therefore respects chiral symmetry which leads to a considerable 
simpliflcation. 

At the scale kq> the propagator G and the wave function g{—q, q ~ p) should be optimized 
for a most complete elimination of terms quartic in the quark flelds. In the present context we 
will, however, neglect the momentum dependence of Z^^^, Z^ ^ and hk. The mass term fff' was 
found in [31| for the simple truncation (24) with = 1, m = m = to be Tn{k^) ~ 120 MeV. 
In view of the possible large truncation errors we will take this only as an order of magnitude 
estimate. Below we will consider the range m(/c$) = (45 — 120) MeV for which chiral symmetry 
breaking can be obtained in a two flavor model. Furthermore, we will assume, as usually done 
in large- Ai'c computations within the NJL-model, that Z$(A;$) = Z$ ^^,(g = 0) <^ 1. Moreover, 
the quark wave function renormalization Z^{k) = Z^^k{q = 0) is set to one at the scale for 
convenience. For k < k<s, we will therefore study an effective action for quarks and mesons in 
the truncation 

Ffc = 1 d'^xi^Z^'^j^r + tr [9^<Dta^$] + [/fc($, <l>t) 

+ (^^$.5 - ^^($%) (35) 
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with compositeness conditions 



Uk^{^,^^) = m^(A;$)tr<l>^$ - ^z7(A;$) (det^ + det^'^ 

+ -Ai(tr<l>'^$) H ^ — A2tr f <l>^$ - — tr<l>"^$j +... 

m2(fc$) = ~ (45 MeV)2 - (120 MeV)^ 

h{ki^) = Z^{k<s,) = 1 
Z^{k^) < 1 . (36) 

As a consequence, the initial value of the renormalized Yukawa coupling which is given by 
h{k<i>) = h{k^)Z^^{kis>)Z^^^'^{kis,) is large! Note that we have included in the potential an 
explicit Ua{^) breaking term |3 ~ which mimics the effect of the chiral anomaly of QCD to 
leading order in an expansion of the effective potential in powers of $. Because of the infrared 
stability discussed in section |^ the precise form of the potential, i.e. the values of the quartic 
couplings Ai(/c$) and so on, will turn out to be unimportant. 

We have refrained here for simplicity from considering four quark operators with vector and 
pseudo-vector spin structure. Their inclusion is straightforward and would lead to vector and 
pseudo-vector mesons in the effective action (|35D. We will concentrate first on two flavors and 
consider only the two limiting cases V = and u —>■ oo. We also omit first the explicit quark 
masses and study the chiral limit j = 0. (The more realistic three flavor situation with a finite 
but non- vanishing V and non-vanishing quark masses is depicted in section |^.) Because of the 
positive mass term m'^{k<^) one has at the scale a vanishing expectation value ($) = (for 
J = 0). There is no spontaneous chiral symmetry breaking at the compositeness scale. This 
means that the mesonic bound states at k^ and somewhat below are not directly connected to 
chiral symmetry breaking. 

The question remains how chiral symmetry is broken. We will try to answer this question 
by following the evolution of the effective potential Uk from kq, to lower scales using the exact 
renormalization group method outlined in section |^ with the compositeness conditions ( ^6|) 
defining the initial values. In this context it is important that the formalism for composite 
fields also induces an infrared cutoff in the meson propagator. The flow equations are 



therefore exactly of the form (|TT]) (except for the supertrace), with quarks and mesons treated 
on an equal footing. In fact, one would expect that the large renormalized Yukawa coupling 
will rapidly drive the scalar mass term to negative values as the IR cutoff k is lowered |TD[ . This 



will then finally lead to a potential minimum away from the origin at some scale k^sB < kq, 
such that ($) 7^ 0. The ultimate goal of such a procedure, besides from establishing the onset 
of chiral symmetry breaking, would be to extract phenomenological quantities, like or meson 
masses, which can be computed in a straightforward manner from in the IR limit k ^ 0. 

At first sight, a reliable computation of Ffc^o seems a very difficult task. Without a trunca- 
tion Ffc is described by an infinite number of parameters (couplings, wave function renormal- 
izations, etc.) as can be seen if Tk is expanded in powers of fields and derivatives. For instance, 

^^The anomaly term in the fermionic effective average action lias been computed in pTl. 
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the pseudoscalar and scalar meson masses are obtained as the poles of the exact propagator, 
lim^^o rfc^''((3') I #=($), which receives formally contributions from terms in with arbitrarily 
high powers of derivatives and the expectation value (Jq. Realistic nonperturbative truncations 
of Ffc which reduce the problem to a manageable size are crucial. The remainder of this work 
will be devoted to demonstrate that there may be a twofold solution to this problem: 

• Due to an IR fixed point structure of the flow equations in the symmetric regime, i.e. for 
kxSB < k < kq,, the values of many parameters of Tk for A; — > will be approximately 
independent of their initial values at the compositeness scale For large enough h{k<^) 
only a few relevant parameters (m^(/c<i,), v{ki^)) need to be computed accurately from 
QCD. They can alternatively be determined from phenomenology. 

• Physical observables like meson masses, decay constants, etc., can be expanded in powers 
of the quark masses. This is similar to the way it is usually done in chiral perturbation 
theory To a given finite order of this expansion only a finite number of terms of a 



simultaneous expansion of F^ in powers of derivatives and $ are required if the expansion 
point is chosen properly. 

In combination, these two results open the possibility for a perhaps unexpected degree of 
predictive power within the linear meson model. 

We wish to stress, though, that a perturbative treatment of the model at hand, e.g., using 
perturbative RG techniques, cannot be expected to yield reliable results. The renormalized 
Yukawa coupling is expected to be large at the scale and even the IR value of h is still 
relatively big 

/i(A; = 0) = ^ ~ 7.5 (37) 

JIT 

and h increases with increasing k. The dynamics of the linear meson model is therefore clearly 
nonperturbative for all scales k <kq,. 



5 Flow equations for the linear meson model 

We will next turn to the ERGE analysis of the linear meson model which was introduced in 
the last section. To be specific, we will try to attack the problem at hand by truncating F^ 
in such a way that it contains all perturbatively relevant and marginal operators, i.e. those 
with canonical dimensions (ic < 4 in four space-time dimensions and ignore the evolution and 
effects coming from operators with dc > 4. The effective potential Uk is a function of only four 
SUl{N) X SUr{N) invariants ioi N = 3: 

p = tr<l>^<l> 

^^The invariant T3 is only independent for A'^ > 3. For TV = 2 it can be eliminated by a suitable combination 
of T2 and p. The additional ?7/i(l) breaking invariant oj — i{det $ — det is CV violating and may therefore 
appear only quadratically in Uk- It is straightforward to see that is expressible in terms of the invariants 
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r. = ^tr($t$_i_p 

^ = det $ + det . (38) 

We will expand Uk as a function of these invariants around its minimum in the absence of 
explicit SU (3) breaking due to current quark masses, i.e. p = po = Na^, ^ = = and 
T2 = T3 = where 

ao^^tr($). (39) 

Restricting Uk to operators of canonical dimension < 4: therefore this yields in the chirally 
symmetric regime, i.e., for k^^B < k < where = 

Uk = nfp + -Aip2 + -^\2T2 - -Vi (40) 
whereas in the SSB regime for k < k^^B we have 



1^ 



2 



4 -V / - 2 



Uk = i:Xiik) p - Nam ' + ^—r^Hk)r2 + ^-V{k) a^~\k)p - i . (41) 



Before continuing to the actual beta functions for the couplings or parameters still kept 
in Tk it is worthwhile to pause here and emphasize that naively (perturbatively) irrelevant 
operators can by no means always be neglected. The most prominent example for this is QCD 
itself. It is the very assumption of our treatment of chiral symmetry breaking (substantiated 
by the results of 0]) that the momentum dependence of the coupling constants of some six- 



dimensional quark operators (gg)^ develop poles in the s-channel indicating the formation of 
mesonic bound states. On the other hand, it is quite natural to assume that <I>^ or <1>^ operators 
are not really necessary to understand the properties of the potential in a neighborhood around 
its minimum. Yet, truncating higher dimensional operators does not imply the assumption that 
the corresponding coupling constants are small. In fact, this could only be expected as long as 
the relevant and marginal couplings are small as well. What is required, though, is that their 
influence on the evolution of those couplings kept in the truncation, for instance, the set of 
equations ( ^3)) below, is small. In this context it is perhaps also interesting to note that the 
truncation (^) includes the known one-loop beta functions of a small coupling expansion as 
well as the leading order result of the large-iV^ expansion of the Ul{N) x Uji{N) model p9[ . 



This should provide at least some minimal control over this truncation, even though we hope 
that our results are significantly more accurate. 

Inserting the truncation (0), (j^, (0) into (|lT]) reduces this functional differential equation 
for infinitely many variables to a finite set of ordinary differential equations. This yields, in 
particular, the beta functions for the couplings Ai, A2, V and nf' or Wq. Details of the calculation 
can be found in [^]. We will refrain here from presenting the full set of flow equations but 
rather illustrate the main results with a few examples. Defining dimensionless renormalized 



18 



mass, VEV and coupling constants 



e{k) = k-^m\k) = m\k)Z^\k)k~^ 

K{k) = k^-'^Nal = Z^e-^Nwl 

h\k) = h\k)Z^\k)Z:p\k) 

\{k) = XiZ^^{k) ; i = l,2 

N 

v{k) = i7{k)Z^''{k)k^-^ . (42) 
one finds, e.g., for tlie spontaneous symmetry breaking (SSB) regime and 17 = 

1 



TTT = -(2+77$)/t 



'dt 



dt 



+ (iV^-l) 



1 + 



167r2 

A2 

Ai 



167r2 



N'Xiq{0) + 9Xilt{2XiK 



+ {N' - 1) [Ai + X2]'ltiX2^) - A^hHr"' ^ 



N ' 



(43) 



Y^{^Ap^(0) + ^(iV^-4)Ap^(A2«:) 



2r7$A2 

iiV2A^/t_i(0, A2/«) + 3[A2 + 4Ai]A2/ti(2Ai«, A2K) 



N 



{N^ - l)lff^'^{^Kh^,e + KX2] 11^,7]^) -/S5^^'^(^K/i^e + 2KAi; 77^,77,^) I 



Here r]q, = —dt\iaZ<s,, rj^ = —dtlnZ^ are tlie meson and quark anomalous dimensions, respec- 
tively The symbols 4^"*^ denote mass threshold functions. A typical example 
is 



Z„(w) = SriTT k 



2 7 2n-4 



A 9t(Z^^i?,(g2)) 



(44) 



(27r)4 [P(g2) + 

with P{q^) = + Z'^^ Rk{q^) . These functions decrease monotonically with their arguments 
w and decay ~ yj-C'^+i) for w ^ 1. Since the arguments w are generally the (dimensionless) 
squared masses of the model, the main effect of the threshold functions is to cut off quantum 
fluctuations of particles with masses ^ k"^. Once the scale k is changed below a certain mass 
threshold, the corresponding particle no longer contributes to the evolution of the couplings 
and decouples smoothly. 
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Within our truncation the beta functions ( ^31) for the dimensionless couphngs look almost 
the same as in one-loop perturbation theory. There are, however, two major new ingredients 
which are crucial for our approach: First, there is a new equation for the running of the mass 
term in the symmetric regime or for the running of the potential minimum in the regime 
with spontaneous symmetry breaking. This equation is related to the quadratic divergence 
of the mass term in perturbation theory and does not appear in the Callan-Symanzik 



or 



Coleman-Weinberg g2| treatment of the renormalization group. Obviously, this equation is 
the key for a study of the onset of spontaneous chiral symmetry breaking as k is lowered from 
/c$ to zero. Second, the most important nonperturbative ingredient in the flow equations for 
the dimensionless Yukawa and scalar couplings is the appearance of effective mass threshold 
functions like ( ^If ) which account for the decoupling of modes with masses larger than k. Their 
form is different for the symmetric regime (massless fermions, massive scalars) or the regime 
with spontaneous symmetry breaking (massive fermions, massless Goldstone bosons). Without 
the inclusion of the threshold effects the running of the couplings would never stop and no 
sensible limit k ^ could be obtained due to unphysical infrared divergences. The threshold 
functions are not arbitrary but have to be computed carefully. The mass terms appearing 
in these functions involve the dimensionless couplings. Expanding the threshold functions in 
powers of the mass terms (or the dimensionless couplings) makes their non-perturbative content 
immediately visible. (It is the presence of threshold functions which explains why one-loop type 
formulae could be used for a necessarily nonperturbative computation of critical exponents in 
three-dimensional scalar theories lEBl.) 



6 The chiral anomaly and the 0(4)— model 

We have seen how the mass threshold functions in the flow equations describe the decoupling of 
heavy modes from the evolution of Tk as the IR cutoff k is lowered. In the chiral limit with two 
massless quark flavors {N = 2) the pions are the massless Goldstone bosons whereas all other 
mesons have masses larger than and are therefore practically decoupled already at the scale 
The effect of the physical pion mass of ~ 140 MeV, or equivalently of the two small but 
non-vanishing current quark masses, can then easily be mimicked by stopping the flow of at 
k = by hand. This situation changes significantly once the strange quark is included. Now 
the 7] and the four K mesons appear as additional massless Goldstone modes in the spectrum. 
They would artificially drive the running of F^ at scales ^ k ^ 500 MeV where they should 
already be decoupled because of their physical masses. It is therefore advisable to focus on the 
two flavor case = 2 as long as the chiral limit of vanishing current quark masses is considered. 

It is straightforward to obtain an estimate of the (renormalized) coupling u parameterizing 
the explicit ?7a(1) breaking due to the chiral anomaly. From (^Tj) we find 

mj, = ^Jy^o''^ ^ 1 GeV (45) 

which translates for N = 2 into 

i^{k ^ 0) ~ 1 GeV . (46) 
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This suggests that z/ — oo can be considered as a reahstic hmit. An important simphfication 
occurs for = 2 and u ^ oc, related to the fact that for = 2 the chiral group SUl{2) x 
SUr{2) is (locally) isomorphic to 0(4). Thus, the complex (2, 2) representation $ of SUl{2) x 
SUr(2) may be decomposed into two vector representations, (o", vr*^) and {ri',a^) of 0(4): 

$ = i (a - ir]') + ^ (a^ + ^t^^) ■ (47) 

For V oo the masses of the rj' and the diverge and these particles decouple. We are then 
left with the original 0(4) symmetric linear cr-model of Gell-Mann and Levy |^ coupled to 



quarks. The flow equations of this model have been derived previously |]T0|, for the truncation 
of the effective action used here. For mere comparison we also consider the opposite limit i/ — > 0. 
Here the i]' meson becomes an additional Goldstone boson in the chiral limit which suffers from 
the same problem as the K and the rj in the case A^ = 3. Hence, we may compare the results 
for two different approximate limits of the effects of the chiral anomaly: 

• the 0(4) model corresponding to A^ = 2 and v ^ oo 

• the Ul{2) X Ur{2) model corresponding to A^ = 2 and z/ = 0. 

For the reasons given above we expect the first situation to be closer to reality. In this case we 
may imagine that the fluctuations of kaons, rj rj' and the scalar mesons (as well as vector and 
pseudovector mesons) have been integrated out in order to obtain the initial values of — 
in close analogy to the integration of the gluons for the effective quark action rfcp['?/^] discussed 
in section ^. We will keep the initial values of the couphngs m^(fc$), Ai(/c$) and Z^{k^) as 
free parameters. Our results should be quantitatively accurate to the extent to which the local 
polynomial truncation is a good approximation. 



7 Infrared stability 

Eqs. (^) and the corresponding set of flow equations for the symmetric regime constitute a 
coupled system of ordinary differential equations which can be integrated numerically. Similar 
equations obtain for the 0(4) model where N = 2 and the coupling A2 is absent. The most 
important result is that chiral symmetry breaking indeed occurs for a wide range of initial 
values of the parameters including the presumably realistic case of large renormalized Yukawa 
coupling and a bare mass m(/c$) of order 100 MeV. A typical evolution of the renormalized 
mass m{k) is plotted in figure |I|. Driven by the strong Yukawa coupling, m decreases rapidly 
and goes through zero at a scale fc^sB not far below kq,. Here the system enters the SSB regime 
and a non-vanishing (renormalized) VEV ct/j for the meson field $ develops. The evolution 
of (Tiiik) turns out to be reasonably stable already before scales k ^ m.,; where the evolution 
is stopped. We take this result as an indication that our truncation of the effective action 
leads at least qualitatively to a satisfactory description of chiral symmetry breaking. The reason 
for the relative stability of the IR behavior of the VEV (and all other couplings) is that the 
quarks acquire a constituent mass Mq = haR ~ 350 MeV in the SSB regime. As a consequence 
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Figure 1: Evolution of the renormalized mass m in the symmetric regime (dashed hne) and 

1 /2 

the vacuum expectation value aR = ao of the scalar field in the SSB regime (solid line) as 
functions of k for the Ul{2) x Ur{2) model. Initial values are Aj(fc$) = 0, K^iki^) = 300 and 
m{k^) = 63 MeV. 

they decouple once k becomes smaller than Mg and the evolution is then dominantly driven 
by the massless Goldstone bosons. This is also important in view of potential confinement 
effects expected to become important for the quarks for k around Aqcd — 200 MeV. Since 
confinement is not explicitly included in our truncation of Tk, one might be worried that such 
effects could spoil our results completely. Yet, as discussed in some more detail in section 
^, only the colored quarks should feel confinement and they are no longer important for the 
evolution of the meson couplings for k around 200 MeV. One might therefore hope that a 
precise treatment of confinement is not crucial for this approach to chiral symmetry breaking. 

Most importantly, one finds that the system of flow equations exhibits an IR fixed point 
in the symmetric phase. As already pointed out one expects to be rather small at the 
compositeness scale In turn, one may assume that, at least for the initial range of running 
in the symmetric regime the mass parameter e ~ Z^^ is large. This means, in particular, that 
all threshold functions with arguments ~ e may be neglected in this regime. As a consequence, 
the flow equations simplify considerably. We find, for instance, for the Ul{2) x Ur{2) model 



dtXi 
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87r2 

This system possesses an attractive IR fixed point for the quartic scalar self interactions 

Ai* = \\2* = ^ . (49) 



Furthermore it is exactly soluble Because of the strong Yukawa coupling the quartic 

couplings Ai and A2 generally approach their fixed point values rapidly, long before the systems 
enters the broken phase (e —>■ 0) and the approximation of large e breaks down. In addition, 
for large the value of the Yukawa coupling at the scale kg where e vanishes (or becomes 

small) only depends on the initial value 

= T^TTTT = • (50) 



h\k^) k 



Hence, the system is approximately independent in the IR upon the initial values of Ai, A2 
and /i^, the only "relevant" parameter being Q e^o- In other words, the effective action looses 
almost all its "memory" in the far IR of where in the UV it came from. This feature of the 
flow equations leads to a perhaps surprising degree of predictive power. In addition, also the 
dependence of f^, = 2cr/j(fc — > 0) on eg is not very strong for a large range in eo, as shown 
in flgure |^. The relevant parameter eo can be flxed by using the constituent quark mass 
Mq = {haR){k = 0) ~ 350 MeV as a phenomenological input. One obtains for the 0(4) model 

eo ~ 0.02 . (51) 

The resulting value for the decay constant is 

= (91 - 100) MeV (52) 

for h?{kq,) = 10^ — 300. It is striking that this comes close to the real value /t^ = 92.4 MeV 
but we expect that the uncertainty in the determination of the compositeness scale ki^, and the 
truncation errors exceed the influence of the variation of /i^(A;$). We have furthermore used 
this result for an estimate of the chiral condensate: 

i^^) = -ioUklZ^'/^k = 0)A ^ -(195 MeV)3 (53) 



where the factor A ~ 1.7 accounts for the change of the normalization scale of (tpip^ from /c<j, 
to the commonly used value 1 GeV. Our value is in reasonable agreement with results from 



^^Once quark masses and a proper treatment of the chiral anomaly are included for N — 3 one expects that 
rriq and v are additional relevant parameter. Their values may be fixed by using the masses of tt, K and 77' as 
phenomenological input. 
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Figure 2: The pion decay constant as a function of eo for kq> = 630 MeV, Xi{k^) = A2 (/::$) = 
and h?{ki^) = 300 (solid line) as well as h?{kq>) = 10"^ (dashed line). 



sum rules |^^. This result is non-trivial since not only /c<i> and enter but also eo and the IR 
value Z$(0). Integrating (^) for 77$ one finds 

Z^{k) = Z^{k^) + ^In'^ . (54) 

Thus Zi^{k) will indeed be practically independent of its initial value Z<i,(fc$) already after some 
running as long as Zq,{kq,) is small compared to 0.01. 

The alert reader may have noticed that the beta functions (|48| ) correspond exactly to those 
obtained in the one-quark-loop approximation or, in other words, to the leading order in 



the large-iVc expansion for the Nambu-Jona-Lasinio model The fixed point (^) is then 
nothing but the large-A^^^ boundary condition on the evolution of Ai and A2 in this model. Yet, 
we wish to stress that nowhere we have made the assumption that N^. is a large number. On the 
contrary, the physical value Nc = ?> suggests that the large- A^'c expansion should a priori only be 
trusted on a quantitatively rather crude level. The reason why we expect (|48|) to nevertheless 
give rather reliable results is based on the fact that for small Zq,{kq,) all (renormalized) meson 
masses are much larger than the scale k for the initial part of the running. This implies that the 
mesons are effectively decoupled and their contribution to the beta functions is negligible leading 
to the one-quark-loop approximation. Yet, already after some relatively short period of running 
the renormalized meson masses approach zero and our approximation of neglecting mesonic 
threshold functions breaks down. Hence, the one-quark-loop approximation is reasonable only 
for scales close to /c$ but is bound to fail around /c-^sb and in the SSB regime. What is important 
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in our context is not the numerical value of the partial fixed points (^) but rather their mere 
existence and the presence of a large coupling h"^ driving the Aj fast towards them. This is 
enough for the IR values of all couplings to become almost independent of the initial values 
Aj(/c$). Similar features of IR stability are expected if the truncation is enlarged, for instance, 
to a more general form of the effective potential f/fc<j,. 



8 Phenomenology of the hnear meson model for three 
flavors 

In the preceding sections we have presented an investigation of the linear meson model subject 
to some simplifications and approximations: 

• We have restricted ourselves to the case of two light quark fiavors, N = 2, whereas it is 
known from chiral perturbation theory that the strange quark can also be considered 



as light. 

• Explicit chiral symmetry breaking due to current quark masses has been neglected. For 

= 2 its main effect is to give the pions a non-vanishing mass which is easily mimicked 
by stopping the RG evolution around k = by hand. For = 3 this would be rather 
inaccurate, since the kaons and the rj have masses quite different from the pions and 
should therefore decouple much earlier. This necessitates an explicit consideration of the 
source term ~ j. 

• We have treated the effects of the chiral anomaly only for two limiting cases: u = and 
u —>■ oo. In reality u is, however, large but finite. 

• The IR-cut-off effective action F^ was truncated in such a way that only operators of 
canonical dimension dc < 4: were kept. 

It is, of course, necessary to investigate the effects of these simplifications systematically. It is 
straightforward to relax the first three points, i.e. to include three massive quark fiavors and 
allow for a non-vanishing and finite u. A systematic study of the truncation errors is more 
involved. Even though it is possible to include additional terms in F^ it would be desirable 
to have a guiding principle at hand, which would allow to determine the invariants which one 
should keep in a truncation of F^ and those which can be neglected. Clearly, such a principle 
will depend on the kind of physical questions one poses. In the following we will apply three 
different criteria for this purpose: 

1. Phenomenological observables like meson masses or decay constants can be expanded 
systematically in powers of (current) quark masses in the linear meson model in a similar 
way as it is usually done in chiral perturbation theory. This can be used to determine the 
operators of F^ which are necessary to compute an observable to a given order in ruq. 
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2. We have assumed so far that all other low-energy degrees of freedom of strong inter- 
action physics like vector and axial-vector mesons have been integrated out. Carrying 
out this procedure explicitly in some cases induces relatively large values for several 
non-renormalizable interactions in the linear meson model. They have to be taken into 
account, at least for comparison with phenomenology. 

3. Even though the range of evolution from to A; = is not too large {mT^/k^ ^ 0.2), 
it is likely that higher dimensional operators are dimensionally suppressed at A; = m,r as 
long as their anomalous dimension is not very large. This is related to the "triviality" 
of the linear quark meson system, i.e. the existence of the interaction free "Gaussian" 
fixed point. We will therefore assume that of two operators contributing to a physical 
observable at the same order in the quark mass expansion, the lower dimensional one 
dominates, provided the other one does not receive large corrections from other degrees 
of freedom assumed to be integrated out. This can also be understood as an expansion 
in powers of the chiral condensate ~ Wq. 

We will sketch in the following how these criteria can be applied to the linear meson model. 
The details of this procedure can be found in |H5 . 



In order to perform a quark mass expansion of meson observables it is useful to decompose 
$ with respect to the vector-like SUv{S) symmetry 

" ^° 71 i'^" ^ Ts^" ^ ^ 7!^') ■ ^^^^ 

For = 3 the hermitean traceless 3x3 matrices $p and $s represent (up to mixing effects) the 
pseudoscalar and scalar octets, respectively. The real fields Xp and Xs are associated with the 
parity odd and even singlets corresponding to the r]' meson and the cr-resonance, respectively. 
The constant Wq is chosen as the VEV of $ in the limit of equal quark masses rUq = {rriu + 
TTT'd + "T's)/3. Differences between the strange, down and up, quark masses imply 

i^s) = ^ (^A3 - VSvXs) . (56) 

Here A3 and Ag denote the two diagonal Gell-Mann matrices. The parameters w and v corre- 
spond to the explicit isospin and SUv{S) violation in the pseudoscalar meson decay constants 



W = fK± - f 



K° 



1 



^^j ^ ^[fK± + fKo-2f^}. (57) 

Here and denote the wave functions renormalization constants of $p and $s, respectively. 
The /j are related to the physical decay constants fi by additional wave function renormalization 
effects which will be discussed below. 
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It is convenient to expand the effective potential around $ = ctq or p = Sctq, ^ = and 
r2 = Ta = 0: 

f^(p, 7-2,^3,0 = ml{p- po) -^71^ - ^o-ao{p- po)] 

1- 2 1- 1- 

+ (P - Po) + 2-^2^2 + -AsTg (58) 

+ ^^i(p-Po)(e-eo) + ^^2(p-Po)r2 
+ ^^3(e-eo)r2 + ^^4(e-eo)' + ... . 

The bare meson mass squared matrix is obtained from the second derivatives of U with respect 
to the components of the fields $p, Xs and Xp evaluated at the minimum of f/^ — | tr 
Decomposing the chiral invariants p, T2, and C, with respect to the SUvi'i) multiplets 
$p, Xs and Xp and taking into account that the only non-vanishing VEV is ~ 0{mq) 
it is straightforward to see that the potential ( [58| ) contains all terms which are necessary and 
sufficient for a computation of the bare pseudoscalar meson masses to second order and the 
scalar octet and pseudoscalar singlet masses to linear order in the quark masses Q. 

In a similar way one can see that modifications of the minimal kinetic term for $ beyond 
the approximation ( P5| ) are necessary for a determination of the physical meson masses to the 
desired order in m^. There are various two-derivative terms all contributing to the same order 
in the quark mass expansion to, e.g., meson masses. As argued above we will, however, assume 
that the lowest dimensional ones dominate. The ones relevant here are given by 

- -X^\ Tr ($^9.$ - a$"^<l>) _ ^A'^t,!,) 

+ Tr ($9^$^" - d^m^) {^d''^^ - S'^^^t) } 

- ^X^[ Tr ($^9/,$ + d^'i>^<^) (^<^^d^<^ + d^<i>^<^) 

+ Tr {^d^<^^ + 9^$$^) ($9'^$"^ + d^<l><l>^) } + .... 

These modifications of the kinetic term induce different wave function renormalization constants 
Zh, Zm, Zg and Zp for the four multiplets ^s, ^p, Xs and Xp, respectively, e.g. 



Zm 





—2 



Zp — Z^ — 2U^aQ + X^(Tq 
Zh = Z^- U^Wq - X^wl . (60) 



^"^If one counts the effects of tlie chiral anomaly as being of 0{mq), eq. ( [5q ) is also sufficient for the bare 
pseudoscalar singlet or (up to mixing) the rj' meson mass to quadratic order in mq. 
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Furthermore, they lead to a spht in the wave function renormahzations for the individual 
members of each multiplet: 



'Z'-K 

Zk° 



i - u^v 

1 

1 + 2^"i [v + w) 
1 

1 + -t^m [V-W) 
1 + UmV 



where 



-1/2^-1 



The physical meson masses and decay constants are then given by 



Mi 



MlZ-^Z-y, f^ = f^Zl/'; etc. 



(61) 



(62) 



(63) 



where the Mf are the bare meson mass terms given by the eigenvalues of the matrix of second 
derivatives of the potential U. For the pseudoscalar octet the relations ( |63D with (^) contain 
all contributions to the meson masses to second order and the decay constants to linear order 
in the quark mass expansion. 

In addition, the non-minimal kinetic terms ( ^91) result in a correction relevant for the rj-r]' 
mixing. This is related to a quark mass dependence of kinetic terms and off-diagonal kinetic 
terms arising from 



AC 



r^^mZfriZ]^ {Z, ZmZh 



,1/2 



with 



UJ. 



pm 



(2X^cro + U^^ {ZpZhZm) 



(64) 



(65) 



The parameter ujprn induces a momentum dependence in the rj-rj' mixing angle. Consequently, 
the mixing angle will be different for the on-shell decays 77 — 27 and 77' 27 through which 
it is usually determined experimentally. This qualitative result is, in fact, in agreement with 
observation. For a quantitative comparison with experiment it is, however, more convenient to 
compute the measured ratios 



M3 r(7r 



27) 



M3 T{r] 27) 



In 



27) 



M3 r(v ^ 27) 



(66) 



from (1^) and (|^). The results, as well as and M^/, are now functions of f/$ and 
or Zp/Zm and ujm- Using the observed values for M^, Ki^±, Mk", Mj^i, and fx± as input 
we have plotted including all corrections of second order in the quark mass expansion as 
a function of utmV for various values of Zp/Zm in figure ^. Figure ^ shows the ratios (^Up of 
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Figure 3: The plot shows as a function of LOmV for various values of Zp/Z„i. The solid line corresponds to 
Zpl Z„i — 1 and the difference in Zp/Zm between two adjacent lines is 0.1. The horizontal dotted line indicates 
the experimental value Mjj ~ 547.5 MeV. 



decay constants as functions of UmV for the same values of Zp/Zm- The agreement with the 
experimental values is satisfactory for Zp/Zm — 0.9 and UmV — —0.19. In addition, it is clear 
from figure || that a phenomenologically acceptable value of /^/ can only be obtained once 
the non-minimal kinetic terms in ( ^91) are taken into account. We furthermore note that the 
deviation of Zp/Z^ from one is rather small. This is compatible with the assumption that the 
chiral anomaly which induces a non-vanishing f/^ is an effect which is comparable in size to 
corrections linear in the quark masses. 

Interestingly, the non-minimal kinetic term ~ in (|59| ) which is responsible for the 
relatively large value of UmV receives a large contribution from the exchange of axial-vector 
mesons. More precisely, the longitudinal component of i9^p^, with the axial-vector meson 
nonet, represents a state which therefore mixes with the 1] and r]' mesons — a mechanism 
usually referred to as "partial Higgs effect". A rough estimate of the contribution to 
resulting from axial-vector meson exchange yields a value of uJ^^v ~ —0.16. This is already 
quite close to the value UmV ~ —0.19 suggested by figures |] and One is therefore tempted 
to even assume that a// kinetic terms (including those of 0{d^)) beyond the minimal one ~ 
are dominated by the contributions from axial-vector meson exchange. This "leading mixing 
approximation" leads, in fact, to a satisfactory description of the t]--)]' system P5|. 

Analogously one may apply the quark mass expansion to the scalar sector of the linear 
meson model. As an interesting result one obtains a Gell-Mann-Okubo type relation for the 
scalar octet 

Ml = \ (4M^. - Ml) (67) 
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Figure 4: The plots show the ratios fri/f-K and f-q'/f-K as functions of LOmV for various vahics of Zp/Zm- The 
sohd hnes correspond to Zp/Z^ = 1 and the difference in Zp/Z^ between two adjacent hnes is 0.1. The 
experimentally allowed windows {la) for both quantities are bounded by the horizontal solid lines. 
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where M^^ denotes the mass of the scalar partner of the rj meson. This relation can be seen to be 
exact to linear order in the quark mass expansion. Using Mk* = 1430 MeV and Ma^ = 983 MeV 
as input one obtains 

M = 1550 MeV . (68) 



Yet, it can be seen that the relation (pTQ is subject to large corrections of second order in the 
quark mass expansion induced by mixing of the charged scalar states with the longitudinal 
component (9^Py vector meson nonet. 

Finally, it should be mentioned that the linear meson model as outlined so far also leads to 



satisfactory results for the ratios of light quark masses [46|. In particular, one finds 



niu + rus Ml± fK± fZK± 



ruu + ms M|-± fK± (Zk± 



(69) 



Inserting the Zi defined in (|6T|) with uOmV as determined from the rj-rj' sector this leads to 
values for the quark mass ratios which are in excellent agreement with the ones obtained in 
chiral perturbation theory. We take this as strong support for our treatment of the linear 
meson model and, in particular, for the inclusion of the non-renormalizable terms in ( |5BD and 
(^9]). Taking these terms into account appears to be crucial for obtaining results for meson 
masses and decay constants as well as current quark ratios in agreement with experiment or 
well established results from other methods. These terms are the main difference between our 
phenomenological treatment of the linear meson model and that of previous works 



From our phenomenological study we conclude that at least the couplings = tr j/ {2N(Tq) 



V, Ai and A2 should be included in the truncation ( p5| ) of the potential. Also the difference 
between Z^ and Z^ (which turns out to be Z^jZ^ ^ 0.4 for A; — 0) needs to be incorporated. 
This is the minimal setting for which the meson masses have qualitatively the correct values 
for — > 0. It seems likely that a reasonable estimate of ao and Zm and therefore the flavor 
average of the decay constants, i^fx + /7r)/3 — ^a^ZU^, becomes already possible with this 
type of truncation. For a more detailed comparison with experiment the coupling should 
also be taken into account. A computation of A3 would be helpful for an understanding of the 
mass splitting in the scalar (O"*""*") sector. 

We conclude that it would be very desirable to gain additional information on the parameters 
of the linear meson model contained in ( |55D and (|5^. This would allow to establish relations 
between different scalar and pseudoscalar meson observables beyond those based on chiral 
symmetry. In particular, in the scalar sector which is not very well understood this would be 
quite interesting. The application of the exact RG to the effective action (|58|), (|59D, in fact, 
offers this prospect due to the infrared fixed point structure established in section |^. If it is 
indeed realized as indicated by the simplified analysis for two flavors outlined in this work, it 
should allow for a determination of many of the couplings and wave function renormalization 
constants introduced in (pHf) and ([59|). Work in this direction is in progress. 
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